Time-reversal-symmetry-broken quantum spin Hall effect 
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Quantum spin Hall (QSH) state of matter is usually considered to be protected by time-reversal 
(TR) symmetry. We investigate the fate of the QSH effect in the presence of the Rashba spin-orbit 
coupling and an exchange field, which break both inversion and TR symmetries. It is found that the 
QSH state characterized by nonzero spin Chern numbers C± = ±1 persists when the TR symmetry 
is broken. A topological phase transition from the TR symmetry-broken QSH phase to a quantum 
anomalous Hall phase occurs at a critical exchange field, where the bulk band gap just closes. It is 
also shown that the transition from the TR-symmetry-broken QSH phase to an ordinary insulator 
state can not happen without closing the band gap. 
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The quantum spin Hall (QSH) effect is a new topolog- 
ically ordered electronic state, which occurs in insulators 
without a magnetic field. [l| A QSH state of matter has 
a bulk energy gap separating the valence and conduction 
bands, and a pair of gapless spin filtered edge states on 
the boundary. The currents carried by the edge states 
are dissipationless due to the protection of time reversal 
(TR) symmetry and immune to nonmagnetic scattering. 
The QSH effect was first predicted in two-dimensional 
(2D) models [IS]. It was experimentally confirmed soon 
after, not in graphene sheets Q but in mercury telluride 
(HgTe) quantum wells 

Graphene hosts an interesting electronic system. Its 
conduction and valence bands meet at two inequivalent 
Dirac points. Kane and Mele proposed that the intrinsic 
spin-orbit coupling (SOC) would open a small band gap 
in the bulk and also establish spin filtered edge states 
that cross inside the band gap, giving rise to the QSH 
effect Q. The gapless edge states in the QSH systems 
persist even when the electron spin s z conservation is de- 
stroyed in the system, e.g., by the Rashba SOC, and are 
robust against weak electron-electron interactions and 
disorder 0,0]. While the SOC strength may be too weak 
in pure graphene system, the Kane and Mele model cap- 
tures the essential physics of a class of insulators with 
nontrivial band topology @,0|- A central issue relating to 
the QSH effect is how to describe the topological nature 
of the systems. A Zi topological index was introduced 
to classify TR invariant systems Q, and a spin Chern 
number was also suggested to characterize the topological 
order [5j]. The spin Chern number was originally intro- 
duced in finite-size systems by imposing spin-dependent 
boundary conditions Recently, based upon the non- 
commutative theory of Chern number Q, Prodan [l(| 
redefined the spin Chern number in the thermodynamic 
limit through band projection without using any bound- 
ary conditions. It has been shown that the Zi invariant 
and spin Chern number yield equivalent description for 
TR invariant systems [HH12 1 . 



The QSH effect is considered to be closely related to 
the TR symmetry that provides a protection for the edge 
states and the Z^ invariant. An open question is whether 
or not we can have QSH-like phase in a system where 
the TR symmetry is broken. Very recently, it was sug- 
gested [13| that the quantum anomalous Hall (QAH) ef- 
fect can be realized in graphene by introducing Rashba 
SOC and an exchange field. In this Letter, we study the 
Kane and Mele model by including an exchange field. We 
calculate the spin Chern number C s analytically, and use 
this integer invariant to distinguish different topological 
phases in the model with breaking TR symmetry. We 
find a TR symmetry-broken QSH phase with C± = ±1, 
indicating that the QSH state could survive, regardless 
of the broken TR symmetry, until the exchange field is 
beyond a critical value, at which the bulk band gap closes 
and reopens, and the system enters a QAH phase with 
C± = 1 (or —1). By further inclusion of an alternating 
sublattice potential, we show that the transition from 
the TR-symmetry-broken QSH phase to an ordinary in- 
sulator state is generally accompanied by closing of the 
band gap. Our conclusion extends the conditions under 
which the topological QSH state of matter can happen, 
and opens the door to magnetic manipulation of the QSH 
effect. 

We begin with the Kane and Mele model defined on a 
2D honeycomb lattice [H, [|| with the Hamiltonian 

H = -t^2 4 c j + "% zZ 4° ■ (dkj x d lk ) Cj 

<ij> * <U> 

+iV R ^2 4^2 • {cr x dij)cj + g^c\a z a. (1) 

<ij> i 

Here, the first term is the usual nearest neighbor hopping 
term with cj = (cL>cji) as the electron creation opera- 
tor on site i and the angular bracket in standing 
for nearest-neighboring sites. The second term is the in- 
trinsic SOC with coupling strength V SOl where a are the 
Pauli matrices, i and j are two next nearest neighbor 
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sites, k is their unique common nearest neighbor, and 
vector dik points from k to i. The third term stands for 
the Rashba SOC with coupling strength Vr, and the last 
term represents a uniform exchange field of strength g. 
For convenience, we will set H, t and the distance between 
next nearest neighbor sites all to be unity. 

We expand Hamiltonian ([1]) in the long-wavelength 
limit at Dirac points K and K' to the linear order in 

kl + kl B. The base 



the relative wave vector k — ^Jk 2 + k 2 

vectors are chosen as {ca-[, ca\, est, cb\}, with A and B 
standing for the two sublattices. We consider first the 
relatively simple case where g = |^. It is straightfor- 
ward to find that for Vr < V s0 , there is a finite energy 
gap, Ae = V3(V so — Vr), which corresponds to a topo- 
logical insulating state exhibiting the QSH effect. For 
Vr > V so , the gap vanishes, and the conductance and 
valence bands cross at Dirac points K and K' . The wave 
functions for the two valence bands near the Dirac point 
K are given by 



<f>i,2(k) = *i,a(*0 
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Here, E 1>2 = k 2 + {V R ±V so f ± Vr) are the cor- 

responding eigenenergies with subscripts 1 and 2 repre- 
senting two bands originating from the sublattice degrees 
of freedom, Fi(k) and F 2 (k) are normalization constants, 
and 9 is the polar angle of k in the reciprocal space. 

For finite g, the obtained analytic expressions for the 
eigenenergies are too long to write here. The energy gap 
Ae between the conduction and valence bands is plotted 
in Fig. la as a function of |<?|/Ko f° r some different values 
of Vr/V so . It is found that, for Vr < V so with increasing 
\g\ from 0, the gap first decreases; as \g\ reaches a crit- 
ical value g Cl the gap closes at k = point; and as \g\ 
further increases, the gap reopens. The critical exchange 
energy g c is determined by the condition of touching the 
conduction and valence bands. For Vr < V so , we have 
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It indicates that g c decreases with increasing Vr/V so . For 
Vr > V so , we have g c — 0, and the band gap always ex- 
ists for finite g. As will be argued below, the insulating 
state for \g\ < g c corresponds to the QSH phase, while 
that for \g\ > g c is also topologically nontrivial with gap- 
less chiral edge states, exhibiting a quantized charge Hall 
conductance. 

The definition of spin Chern number C s relies on a 
smooth decomposition of the occupied valence band into 
two sectors through diagonalization of the electron spin 
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FIG. 1: (a) Normalized energy band gap Ae/V so and (b) 
spectrum gap of a z as functions of |ff|/14o for some different 
values of Vr/V so - 



operator s z — ^a z in the valence band. [lfj| Since s z com- 
mutes with momentum, the decomposition can be done 
for each k separately. To simply show the calculation pro- 
cedure for C s , we first discuss the case of g = 0, where 
the wave functions <pi(k) and tf2(k) for the valence band 
have been given in Eq. (2). By diagonalizing the 2x2 

matrix 



(ip a (k)\a z \ipp(k)) 
two eigenfunctions of a z as 

1 



with a, f3 = 1,2, we obtain 



V>±(fc) 



V2 



(4) 



The minimal spectrum gap between the eigenvalues of 
a z as a function |<?|/V so for different values of Vr/V so , is 
plotted in Fig. lb. The spectrum gap is always nonzero, 
and so we can unambi guo usly calculate the corresponding 



spin Chern numbers [10 . [11 1 



The spin Chern number 
can be defined as a sum over two Dirac points C± 
Ck± + Ck'±, where for the K point (El [Tol 



Ck'±, where for the K point 
1 
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C K ± = ^ / d 2 kQ K± (k), 



(5) 



with Qx±(k) = ie z ■ (Vfct/>±(fc)| x |VfcV±(fc))- C K>± 
can be defined similarly. By using the polar coordi- 
nate system, it is straightforward to obtain Qx±{k) = 
IJ-P K± (fc) with P K ±(k) = T2F 1 (k)F 2 (k). Substitut- 
ing the expression for Qx±(k) into Eq. (5), we derive 
C K (k>)± to be 



C K {k>)± = [Pk{K')±{oo) - Pk(K')±(°)} 



(6) 
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r < V so , numerical calculation yields 
)±(oo) = and Pk(K')(Q) — -Flj as shown 
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FIG. 2: Calculated P K(K ')±(k) for V R /V to = 0.5. The ex- 
change energy is taken to be (a) g = 0, and (b) g/V ao = 1.2. 
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FIG. 3: Phase diagram determined by the Chern numbers in 
the Vr/Vso versus g/V ao plane. The phase diagram in the half 
plane of Vr/V so < is mirror symmetric to Vr/V so > 0, and 
hence not plotted. 



in Fig. 2a. It then follows Ck{k>)± = ±|j and the total 
spin Chern numbers are C± = ±1. Therefore, at g = 0, 
Vr < V so corresponds to a topological QSH insulator, as 
expected. 

Now we consider the finite g case. As has been dis- 
cussed above, the gaps of energy and spin always exist 
for finite g except for \g\ = <? c , so that the spin Chern 
number can be defined in the whole parameter plane ex- 
cept for \g\ — g c . Using a procedure similar to that in 
the g = case outlined above, we obtain Px±(k) and 
P K ,±(k) for both regions \g\ < g c and \g\ > g c . It is 
found that for \g\ < g c , the curves of PK(K')±(k) are 
very similar to that in Fig. 2a with Pk(K')±(°°) = 
and Pk(K')(0) = T 1 unchanged, yielding C± — ±1. Out 
of this region, we obtain C± = 1 for g > g c (see Fig. 
2b), or C± = — 1 for g < —g c . Figure 3 shows a phase 
diagram determined by spin Chern numbers in the g/V so 
versus Vr/V so plane. There are three topologically dis- 
tinct phases characterized by C± = ±1, C+ = C_ = 1, 
and C-t- = C_ — — 1, respectively. From our calculation, 
the boundary between the different topological phases is 
just the condition of closing the band gap. 

To study the edge states in each region, we calculate 
the energy spectrum of a long ribbon with zigzag edges 



and 240 zigzag chains across the ribbon. For \g\ < g c , 
corresponding to C± — ±1 region in the phase diagram, 
the energy spectrum is shown in Fig. 4a. One can easily 
distinguish the edge states from the bulk states. There 
is a small energy gap in the edge modes as can be seen 
from the inset of Fig. 4a, due to the absence of TR and 
inversion symmetries. At a given Fermi level in the band 
gap, there exist four different edge states labeled as A, B, 
C, and D. Through the analysis of the spatial distribution 
of the wave functions, one can find that states A and B 
localize near one boundary of the ribbon, while C and D 
localize near the other boundary. Take states A and B on 
one boundary for example. From the slope of dispersion 
curves at points A and B, it is easy to determine that the 
two edge states are counterpropagating. We also examine 
the spin polarization of the wave functions, state A being 
almost fully spin-up polarized, and state B spin-down 
polarized. Therefore, in the C± — ±1 region there exist 
two counterpropagating edge states with opposite spin 
polarizations on a sample edge, which give rise to no net 
charge transfer but contribute to a net transport of spin. 

The characteristic of the edge states in the C± = ±1 
region with g =/= discussed above is very similar to that 
for the QSH phase at g — protected by the TR sym- 
metry 0. In particular, they have the same spin Chern 
number C± = ±1, indicating that they belong to the 
same topological class. As a result, we call it the TR 
symmetry-broken QSH phase. For the QSH phase pro- 
tected by the TR symmetry, nonmagnetic impurities do 
not cause backscattering on each boundary, and the spin 
transport in the edge states is dissipationless at zero tem- 
perature. In the TR symmetry-broken QSH phase, there 
is usually a weak scattering between forward and back- 
ward movers, as evidenced by the small energy gap in 
the edge state spectrum, leading to a low-dissipation spin 
transport. 

Similar analysis can be applied to the C± — 1 region, 
where the total Chern number of the filled bands sum 
upto C = 2, corresponding to a QAH phase [3] with 
Hall conductivity a xy = 2e 2 /h. The related edge state 
spectrum is shown in Fig. 4b. It is found that states A 
and C localize at one boundary and propagate along the 
same —x direction, while states B and D localize at the 
other boundary and propagate along the same +x direc- 
tion. As a result, in the QAH phase, two edge states at 
each boundary lead to spin-up and spin-down currents 
propagating along the same direction, yielding a quan- 
tized charge conductance. The symmetry-broken QSH 
and QAH phases are topologically distinct. The topolog- 
ical phase transition between them can occur at |<?| = g c 
where the band gap just closes. 

To further investigate the transition from the TR- 
symmetry-broken QSH phase to an ordinary insula- 
tor state, we include an alternating sublattice potential 
M y~l ir Tc\ci into Hamiltonian (1) with r = ±1 for i on 
sublattice A and B, respectively. For g ^ and M ^ 0, 
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FIG. 4: Energy spectrum of a zigzag-edged graphene ribbon. 
The parameters are chosen to be V ao = 0.1, Vr = 0.05, and 
g = 0.03 (a) and g = 0.15 (b). At a given Fermi level in 
the band gap there exist four different edge states, which are 
labeled as A, B, C, and D. 
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FIG. 5: (a) Phase diagram on the g/g c versus M/V so plane 
for g c — 0.5 * ^V ao , and (b) calculated Chern numbers and 
band gap as functions of M/V so for g — 0.5g c and V so = 0.1. 



since both the TR and two-fold rotation symmetries are 
lifted, the two Dirac corns at K and K' become asym- 
metric, leading to an indirect minimal band gap A# be- 
tween the two Dirac points. We find that for a system 
initially in the QSH phase of \g\ < g c with g c given by 
Eq. (3), as \M\ is increased, the indirect band gap closes 
at a smaller critical value \M\ = M_ and reopens at a 



greater critical value \M\ = M + with M T = g c T Iffl- 
The conduction and valence bands overlap in between, 
i.e., A E = for Af_ < \M\ < M + . The spin Chern 
numbers are calculated from the lattice model by pro- 
jecting the two valence bands into two spin sectors, in a 
similar manner to that shown above for the continuum 
model. (For the phase diagram Fig. 3, numerical calcula- 
tions based upon the lattice model have been performed, 
the obtained result being found to agree with that from 
the continuum model.) The spin Chern numbers are well 
defined only in the regions \M\ < M_ and \M\ > M + , 
where A E > 0. We find C± = ±1 for \M\ < M_ and 
C± = for | M | > M_|_. The phase diagram in the g 
versus M plane obtained is plotted in Fig. 5a. The cal- 
culated C± and minimal band gap A g varying along the 
arrowed dashed line in Fig. 5a are shown in Fig. 5b as 
functions of M/V so . A general feature of the phase dia- 
gram is that the transition between the QSH phase with 
C± = ±1 and ordinary insulator state with C± = is al- 
ways companied by closing of the band gap, which serves 
as another signature that the TR-symmetry-broken QSH 
phase is topologically nontrivial and distinct from an or- 
dinary insulator. 
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